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ABSTRACT 

We consider the lifecycle of dust introduced into the hot interstellar medium in iso¬ 
lated elliptical galaxies. Dust grains are ejected into galactic-scale cooling flows in large 
ellipticals by normal mass loss from evolving red giants. Newly introduced dust rapidly 
enters the hot interstellar plasma and is sputtered away by thermal collisions with ions 
during the slow migration toward the galactic center in the cooling flow. Before the grains 
are completely sputtered away, however, they emit prodigious amounts of infrared radi¬ 
ation which may contribute to the large far infrared luminosities observed in ellipticals. 
The infrared emission depends critically on the sputtering rate. Since our understanding 
of both the plasma and radiation environments in ellipticals is quite good, these galaxies 
provide an excellent venue for studying the physical processes of dust grains and perhaps 
also their composition and size distribution. 

In order to study the global properties of grains in ellipticals we construct a new series 
of King-type galactic models which are consistent with the fundamental plane, galactic 
mass to light ratios and other relevant observational correlations. We describe a new 
“continuity” procedure to construct simple time-dependent gas dynamic models for cooling 
flows. 

Although grains can flow a considerable distance from their radius of origin in the hot 
interstellar medium of some galaxies before being sputtered away, we show that the grain 
size distribution at every radius is accurately determined by assuming in situ sputtering 
of dust grains, completely ignoring advection. This occurs since the stellar density prohle 
is so steep that the majority of grains at any galactic radius is produced locally. 

Although thermal sputtering destroys the grains, we show that the dominant source of 
grain heating is absorption of starlight; grain heating by collisions with energetic thermal 
electrons or X-ray absorption are negligible. Previous studies have claimed that the loss of 
thermal energy from a hot, dusty plasma is dominated by grain heating via electron-grain 
collisions and subsequent IR radiation. However, we show that when self-consistent grain 
sputtering is included the dust-to-gas ratio is reduced and radiative cooling, not electron- 
grain interactions dominates plasma cooling, even for the most massive ellipticals. This 
conclusion is insensitive to the grain size distribution assumed for the stellar ejecta. 

1. INTRODUCTION 

Elliptical galaxies, long thought to be simple ensembles of non-interacting stars, are 
now known to be complicated multicomponent systems of stars and gas. The interstellar 
medium (ISM) in ellipticals is dominated by hot gas which is slowly cooling by radiative 
losses; as a result the ISM gradually flows inward as a cooling flow. X-ray observations 
of optically bright ellipticals reveal gas masses of ~ 10® — IO^Mq at high temperatures 
~ lO^K (Forman et al. 1979; Nulsen, Stewart, & Fabian 1984; Forman, Jones, & Tucker 
1985; Canizares, Fabbiano, & Trinchieri 1987). Radio and optical observations (e.g., 21 
cm, Ha, and CO) indicate that elliptical galaxies typically contain 10^ — 10®Mq of cold 
gas usually located in the core regions (see e.g. Kormendy & Djorgovski 1989 for a review). 

In view of the high temperature of the ISM it is remarkable that elliptical galaxies 
also contain substantial amounts of dust. Optical images of ~ 50% of all ellipticals ex- 
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Mbit dust lanes, disks, or patches (e.g., Sadler & Gerhard 1985; Kormendy & Stauffer 
1987; Ebneter et al. 1988). Since dusty disks would be difficult to see when viewed face 
on, it is possible that all ellipticals contain appreciable dust. Further evidence for dust 
are the large far-infrared IRAS luminosities detected in approximately 50% of ellipticals. 
The ratio of flux energy at 100/um to optical B-band emission has a cosmic spread, 
piqqF{IQQ iim)/ peF^ub) = 0.006 — 0.088 (Jura et al. 1987); for cD galaxies this ratio is 
about ten times larger (Bregman, McNamara, & O’Connell 1990). 

The geometrical distribution of dust and the associated infrared emission in ellipticals 
may be either concentrated near the galactic core or distributed throughout the stellar 
system. Concentrated patches or disks of dusty gas observed near the galactic cores could 
arise either (i) from merger events with galaxies having large amounts of cool, dusty gas 
and/or (ii) from cold gas deposited by galactic cooling flows as the hot ISM slowly radiates 
away its thermal energy. These two possibilities may be distinguished observationally. 
Because of the destruction of grains in the hot ISM by sputtering, the dust to gas ratio 
is expected to be much lower, at least initially, in cool gas condensed from cooling flows. 
Moreover, gas deposited by galactic cooling flows should exhibit the same global rotational 
dynamics and sense of rotation as that of the stellar system (Kley & Mathews 1995). 

In addition to the central dust patches, more widely distributed dust is continuously 
generated by stellar mass loss from evolving red giants throughout the galaxy. Knapp, 
Gunn, & Wynn-Williams (1992) have shown that the relatively hot dust that emits 12|Um 
radiation in ellipticals is spatially distributed like the galactic stars and may therefore be 
circumstellar. However, because of the dynamical interaction of the dusty stellar ejecta 
with the hot ISM, the time that recently ejected gas and dust lingers near the parent 
star is limited. When the dusty stellar ejecta interacts with the hot T Wk ism, 
numerous hydrodynamic instabilities are expected (Mathews 1990) resulting in breakup 
and deceleration by drag forces as the ejecta trails off behind the star. The fragmented 
new gas is rapidly heated by thermal conduction as it merges with the hot ISM, even in 
the presence of a magnetic field. It is easy to show (Mathews 1990) that the cool stellar 
ejecta is separated from its parent star and melts into the ISM on time scales 10^ - 10® 
yrs, so fast that little or no grain destruction can occur until the dust is immersed in the 
hot ISM. Once in contact with the hot ISM of number density n, grains with radius a are 
sputtered away by thermal collisions with protons and helium nuclei; the sputtering lifetime 
is tap ~ 10^(a///m)(n/10“^cm“^)“^ yrs provided ISM temperatures are T > 2 x 10® K 
(Tielens et al. 1994; Draine & Salpeter 1979; Seab 1987). However, during their short 
lifetimes as they sputter away and are advected inward by the galactic cooling flow, the 
grains can nevertheless emit a huge amount of far IR radiation (Dwek 1986; 1987). But 
the dust-to-gas ratio and IR emissivity from heated grains is greatly reduced by the short 
sputtering lifetimes so emission and sputtering must be treated in a self-consistent manner. 
Emission from this distributed dust may in fact contribute substantially to the observed 
IRAS luminosities. Unfortunately the spatial resolution of IRAS observations (3' x 5' at 
100//m) is insufficient to discern whether the emission comes from the dusty clouds near 
the galactic centers or is distributed throughout the galaxies. 

In this and a subsequent paper we study the physical properties and emission from 
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the distributed dust component. The infrared emission from distributed dust associated 
with stellar mass loss must be present in every elliptical and is a non-trivial lower limit to 
the observed fluxes. Fortunately the dust enviro nm ent - particularly the gas density and 
temperature in the ISM - is well known from successful hydrodynamic models of cooling 
flows (e.g. Loewenstein & Mathews 1987). The local gas density depends on the rate of 
mass loss from an aging stellar population a^{t)p^{r) g cm“^ s“^ which is well known and 
insensitive to the stellar IMF (Mathews 1989). While grains are destroyed by ion impacts 
from the hot plasma, we hnd that they are heated primarily by the local radiation held. 
The radiation within ellipticals is also very well known in both the optical-IR (Renzini 
& Buzzoni 1986; Bruzual 1985) and in the near UV as observed with the Hopkins UV 
Telescope (e.g. Ferguson & Davidsen 1993). Soft X-ray and far-UV radiation from the 
hot ISM can be determined either theoretically from an atomic emissivity code (Raymond 
1991) or from direct observation (e.g. Trinchieri, Fabbiano, & Canizares 1986; Serlemitsos 
et al. 1993). Our excellent knowledge of the radiation and plasma environment makes 
the interstellar medium in ellipticals an excellent laboratory for studying the physics of 
grain sputtering and emission. In comparing theoretical and observed infrared emission 
the primary uncertainties are the initial distribution of dust grain sizes and the grain 
composition in stellar ejecta; indeed, information about these parameters can be found by 
detailed comparison of theoretical models with infrared observations. A detailed study of 
the expected IR spectrum from ellipticals and their emission into each IRAS band will be 
the subject of a subsequent paper. 

In this paper we discuss the life cycle, thermal properties, and size distribution of 
dust grains ejected into elliptical galaxy cooling hows (ISM) by evolving stars. In order 
to predict the infrared luminosity due to dust emission from a galaxy of given it is 
necessary to know the sputtering rate, the grain heating and cooling rates, the temperature 
and infrared emissivity of grains of all sizes and at all galactic radii. Since we want to 
explore the variation of infrared to optical luminosities among ellipticals of various masses, 
we develop in §2 a scheme for generating models of the radial stellar distribution within 
elliptical galaxies that is consistent with the fundamental plane and other observed galactic 
properties. In §3 we describe a new procedure to develop simple time-dependent models 
for the cooling hydrodynamics in the model galaxies. These cooling hows are used to study 
the ehects of sputtering and advection on grain size distributions in §4. In §5 we compare 
various grain heating mechanisms: by ambient radiation, by electron and ion impacts, 
and by X-ray absorption. We conclude that grain heating by starlight dominates all other 
heating sources. Grain size distributions for each model galaxy are presented in §6. In §7 
we consider the ehects of the presence of grains on the temperature of the ISM; we show 
that cooling of the hot ISM by thermal collisions with grains is generally small compared 
to optically thin radiative losses. 

2. SELF-CONSISTENT MODELS FOR ELLIPTICAL GALAXIES 

For analytic simplicity and for comparison with previous theoretical studies we assume 
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spherical, King-type model galaxies having stellar density prohles given by 


P*(0 = P*o(l i = rlrc 


(1) 


where Tc is the stellar core radins. The stellar distribntion is terminated at an onter 
bonndary rt, or = ft/rc- From eq. (1), the stellar mass within radins ^ in nnits of the 
stellar core mass M*c = is 


Hr 


(0=3 ln[e + (l-Fe")2]-^(l + ^^) 


and the mass snrface density 


SOA) 


2p*orc O 1 + A^V 

(1 + AOV l + 


( 2 ) 


( 3 ) 


depends on the projected radins A = R/vc- 

The parameters that determine the form of the stellar distribntion within each galaxy 
^ P*o, Tc and Tt - mnst be consistent with observational constraints on the properties of 
ellipticals, particnlarly the fnndamental plane relation (e.g. Dressier et al. 1987; Djorgovski 
& Davis 1987). In Appendix A we describe the procednre we have nsed to derive the stellar 
properties of a one parameter family of galaxy models. Relevant properties of the three 
representative galaxies nsed in onr stndy of dnst evolntion are listed in Table 1. Note that 
the listed ellipticals extend from the brightest only to those of intermediate Inminosity 
(Lb < 10^^ Lq) for which the presence cooling flows can be reasonably assnred. The 
global X-ray emission from ellipticals of lower Inminosity is dominated by stellar sonrces, 
evidently low mass X-ray binaries. In addition, low Inminosity ellipticals have shallower 
potential wells and the ISM can become a galactic wind if the collective Inminosity of 
Type la snpernovae exceeds the rate of emission from the hot ISM. We assnme that most 
or all moderate to massive ellipticals have cooling flows as evidenced by their large X-ray 
Inminosities Lx- Recent observations of low iron abnndances in several bright ellipticals 
(Serlemitsos et al. 1993; Forman et al. 1993; Mnshotzky et al. 1995) indicate that the 
Type la snpernova rate is snfhciently low to permit cooling flows in all galaxies considered 
here (Loewenstein & Mathews 1987; 1991). 

In order to estimate the spattering rate of grains in the ISM, it is necessary to have 
estimates of the gas temperatnre, density and velocity in the cooling flows. We describe 
how the gas density and velocity are determined in §3. The temperatnre, ideally, shonld 
be determined by observations of the X-ray emitting gas in galaxies having a spread in 
Inminosities. However, cnrrent temperatnre determinations are available for only a few 
galaxies (e.g., NGC 4472 and NGC 1399: Serlemitsos et al. 1993; NGG 4636: Mnshotzky 
et al. 1995) and are restricted to the highest Inminosity cases. We therefore adopt a 
theoretical approach for determining the temperatnres. Gas ejected by stars in the galaxy 
is rapidly virialized to the mean stellar temperatnre (T*) which characterizes the average 
stellar velocity dispersion (weighted by the stellar mass density) as stars orbit in the total 
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galactic potential. Assuming most of the interstellar gas comes from stars and that other 
sources of heating are small, the temperature of the gas should be very close to (T*). This 
has been conhrmed by detailed hydrodynamic calculations. To estimate (T*) we adopt 
appropriate density structures for the galactic dark matter halos and solve the equation of 
stellar hydrodynamics. This is described in Appendix A. The temperatures listed in Table 
1 are computed in this ma nn er. 

As a check of the above procedure, we plot in Figure 1 the (Lb,(T*))- relation implied 
by our procedure for determining the gas temperature. The location of several of our 
model galaxies are indicated by hlled circles and the general relation is given by the solid 
line. Three galaxies with recent gas temperature determinations by ASCA, NGC 4472, 
NGC 4406, and NGG 4636 (Makishima 1994; Mushotzky et. al. 1995) are also plotted. 
The observed temperatures of these galaxies agree reasonably well with our predicted 
(LB,(T*))-relation. A more accurate determination of the temperatures is not required for 
the considerations of this paper. 


3. TIME-DEPENDENT “GONTINUITY” MODELS FOR GALAGTIG GOOLING FLOWS 
For three nearby elliptical galaxies which were well resolved by the Einstein Observa¬ 
tory Imaging Proportional Gounter (IPG), Trinchieri, Fabbiano, & Ganizares (1986) found 
that the X-ray surface brightness distribution is very similar to that of the starlight, i.e. 
E 2 ;(A) oc S*(A). Since the X-ray emissivity is proportional to the square of the gas density, 
it follows from eq. (1) and the radial constancy of the stellar mass to light ratio that the 
gas density varies as 

p(0 = Po(l + a-"- (4) 

We have implicitly assumed that the gas temperature is independent of galactic radius. 
Isothermality throughout most of the volume of the cooling flow is a general theoretical 
result (Loewenstein & Mathews 1987; 1991) and is supported by recent BBXRT and 
ASCA observations (Serlemitsos et al. 1993; Makishima 1994; Mushotzky et al. 1995). 
Within the galactic cores, however, the gas must cool rapidly and large deviations from 
isothermality are expected. The gas temperature may also decrease at the very largest 
radii since T*(rt) —0 if the stellar velocities become preferentially radial at large radii. 

Since the gas density distribution p(^) of eq. (4) has the same core radius and total 
radius as the stellar distribution, p(^) is completely specihed by the central density po- 
This can be determined by relating the gas density to the X-ray luminosity in the Einstein 
IPG band (~ 0.2 to ~ 4 keV), 


— 1 = 0.165 

nir 


A 


AE 


10 2^ergs cm^ s ^ 


1041 ej-gg g 


-1 


kpc 


(5) 


[ln(2{,)-l]- 


cm 


-3 


where nip is the proton mass. The atomic emissivity in the given energy band, Aae, 
which we have determined from the code of Raymond (1991), depends only on the gas 
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temperature (although this depeudeuce is weak for the temperatures of iuterest). The 
X-ray lumiuosity is observed to correlate with the B-baud optical lumiuosity 

(Douuelly, Faber, & O’Couuell 1990) where h is the Hubble coustaut uormalized to 100 
km s“^ Mpc“^. We assume h = 0.75 throughout. 

Haviug specified the gas deusity distributiou, we ueed to determiue the velocity of the 
gas iu the cooliug flow u{r) to complete the descriptiou of the model galaxies. Although 
the iuward gas flow velocity iu galactic cooliug flows is highly subsouic, these flows are 
uot iu steady state siuce the source term for uew gas expelled by evolviug stars varies 
cousiderably duriug the time required for gas to flow across the galaxy. The specific mass 
loss rate from a siugle-burst populatiou of evolviug stars is a*(t) = an{t/tnow)~^ where 
an ~ 1.6 X 10“^^yr“^, p = 1.35 aud tnow is the preseut time. The time variatiou of a*(t) 
is iuseusitive to the IMF assumed for the galaxy (Mathews 1989). 

Cooliug flows were uot preseut at the very earliest times iu the history of massive 
ellipticals. The siguificaut irou abuudauce iu clusters of galaxies suggests that a cousider- 
able amouut of euriched mass processed iu stars was ejected loug ago from the elliptical 
galaxies (or their progeuitors) by galactic wiuds driveu by Type II superuovae (Mathews 
1989; David, Formau, & Joues 1991). It is therefore reasouable to assume that the cooliug 
flows observed iu massive ellipticals today begau at some time tcf iu the past wheu the 
galactic outflow subsided. 

We uow show that, giveu several simplifyiug assumptious, the velocity field of the 
cooliug flow at the preseut time u{r, tnow) cau be fouud from the coutiuuity equatiou 
aloue, 

% + = (7) 

The iuward flow velocity u{r, t) is the ouly uukuowu depeudeut variable if we uote that the 
stellar deusity p* does uot vary siguificautly with time, aud require that the gas deusity 
at time t vary as 

p{0, ( 8 ) 

where p(^) is the gas deusity at tnow giveu by eq. (4) aud s ~ p/2. Couditiou (8) is 
based both ou the results of detailed hydrodyuamical calculatious (Loeweusteiu & Mathews 
1987; David, Formau, & Joues 1991) aud ou the approximate cousideratious described iu 
Appeudix B. Duriug the evolutiou of cooliug flows the X-ray lumiuosity varies with the 
stellar mass loss rate, oc a*(t). But at auy time is proportioual to the square of the 
gas deusity at auy radius, ~ A as ■ If the gas temperature (aud also Aas) remaius 
roughly uuiform, it follows from these two relatious that p oc a*(t)^/^ oc as iu eq. (8) 
(see Appeudix B). 


7 



Substituting eqs. (1) and ( 8 ) into eq. (7), we get an equation containing only a spatial 
derivative, 

‘ = S-'A.K) (9) 

where 6 = t/tnow: 7 = Po/CinP*otnow is a dimensionless parameter, and the fluid variables 
are expressed in non-dimensional form: 


A. «) 


^ = (i+a‘b 


( 10 ) 


and 


A(0 = ^ = (1 + ^")' 


po 


ri{^.0) = 




Un — 


Ur. 


^ c^nP^o 
po 


( 11 ) 

( 12 ) 


Eq. (9) is valid at all times after the onset of the cooling flow (tc/)- It must therefore 
be valid at the present time, where 6 = 1. On integrating over we determine the velocity 
held at the present time. 


= [A* + syA] (13) 

where no flow across the outer boundary of the galaxy is assumed. 

Eq. (13) represents an approximation for the velocity held in a fully time dependent 
galactic cooling flow. Velocities computed from eq. (13) for the three galaxies of Table 1 are 
shown in Figure 2. The smallest galaxy has the smallest velocity at any given value of the 
normalized radius because the value of 7 is least for the smallest galaxy. The velocity 
profiles compare favorably with those derived from detailed hydro dynamical calculations 
(e.g. Loewenstein & Mathews 1987) as do the flow times predicted from the velocities (see 
Figure 3). 


4. GRAIN SPUTTERING 

Dust grains introduced into hot gas are eroded by collisions with energetic protons 
and helium ions in the gas. The rate at which dust particles are sputtered away has been 
computed in detail by Tielens et al. (1994) and by Draine & Salpeter (1979) for both 
graphite and silicate grains. We adopt the analytic form 


da 

dt 




(14) 


for the rate at which the radius of the dust grain a decreases with time in a hot plasma 
of temperature T and density p. This relation is a good approximation to the detailed 



calculations for both graphite and silicate when h = 3.2 x 10 ^®cm^ s w = 2.5, and 
Td = 2x 10® K. 

Eq. (14) can be used to compute the local sputtering time for a grain, dehned by 


(15) 


tgp ft 


da 

dt 


which can be compared to the flow time to the core radius, dehned by 


tflow (^) fc 



(16) 


These timescales are plotted in Figure 3 for the three galaxy models of Table 1. The 
sputtering time is computed assuming a large grain with a radius of 0.3//m. Smaller 
grains have correspondingly shorter sputtering times. Figure 3 indicates that sputtering 
generally occurs faster than the time for the grains to how to the core in larger galaxies, 
but in smaller galaxies grains may be carried inward by the cooling how for a signihcant 
distance before sputtering away. This suggests that the grain population at a given radius 
in the smaller galaxies will have an admixture of grains that originated further out in the 
galaxy. In contrast, the largest galaxy produces and sputters grains on the spot, except at 
the innermost radii where large grains may how somewhat. 

In the case of the largest galaxy (model a), the grain size distribution in the ISM can 
be determined approximately by balancing the local rate of grain input from stars and the 
rate of grain destruction by sputtering. That is. 


da \ dt J 


S{a) 


(17) 


where N is the number density of grains per unit grain radius, S represents the stellar 
source of grains, and a = da/dt is given by eq. (14). If, for example, we assume that the 
grain size distribution within recently introduced stellar ejecta varies as a power law, i.e. 
S (a) = Aa~^, then eq. (17) is readily solved to give 


N (ft) — N (^amax 




1-9 

max 


-a^-3) 


(18) 


where a 


max 


(14) by 


is the adopted size of the largest grain. The constant h is related to h of eq. 



irip 



(19) 


where T is the (uniform) temperature of the hot ISM. Since we assume that there are no 
sources of grains bigger than aman N{amax) = 0 in eq. (18). Thus, for small grains, the 
size distribution in the ISM goes approximately as a power law in the grain size, N (a) cx 
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but with an exponent that differs by unity from that of the source function. As a 
approaches amaxi the size distribution steepens and falls towards zero. Note that because 
the sputtering rate depends only weakly on T for the temperatures under consideration, 
both d and the grain size distribution are accurately determined despite uncertainties in 
our evaluation of the gas temperature (Figure 1). 

Eq. (18) demonstrates that sputtering significantly modifies the grain size distribution 
from that in the stellar ejecta; smaller grains are relatively less numerous in the intergalactic 
medium compared to that of the stellar ejecta. The general case where inward advection 
of the dust becomes important for determining the grain size distribution will be treated 
in §6, however, the approximate distribution of eq. (18) applies for a large region of galaxy 
model a and the outer regions of galaxy models b and c. The predicted IR spectrum 
and bolometric IR luminosity are greatly modified due to the reduction in the dust-to¬ 
gas ratio and the altered grain size distribution as a result of sputtering. The relative 
underpopulation of smaller grains implies decreased emission at smaller IR wavelengths. 

5. SOURCES FOR DUST HEATING 
5.1 Heating by Ambient Starlight 

Absorption of ambient optical radiation by grains is an important means of grain 
heating. In order to estimate this heating we must compute the mean intensity of starlight 
at any radius in the model galaxies. The mean intensity at radius r is given by 

J,. (r) = ^ y (r, n) = ^ J due (20) 

where ne = cos 9, fi is the solid angle, and is the specific intensity of starlight 
(ergs s“^ cm“^ ster“^ Hz“^). To compute we assume that the stellar emissivity ji, 
(ergs s“^ cm“^ Hz“^) is proportional to the stellar density: 

jAO =3o + (21) 

where jo is an as yet unspecified constant which carries all the frequency dependence of jV- 
We further assume that dust distributed throughout the galactic volume is optically thin to 
stellar light. (This latter assumption can be checked for consistency once the dust content 
of the interstellar gas is computed in §6.) The specific intensity (^, ne) in a direction 
specified by the azimuthal angle ne is computed as an integral along a ray starting at the 
current location and extending to the outer boundary of the galaxy, 

^ ^max 

Iu{i,de)=rc / ju{i) dl. (22) 

Jo 

In this expression, I is a dimensionless length along the line of sight, 

= ^2 ^ ^2 ^ 2^ / cosO, (23) 
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and 


(24) 


^max — + \l ~ 

Eq. (22) can be evalnated to yield 


lu (^, l^e) = jo Tc 


^max + 


(1 + - ei4) A/l + ^" + 2e/f^e/max + /i 

^^0 


{i+e-e^4)V^^ 

The mean intensity of starlight is then given by eq. (20), 


(25) 


J. (0 = 


JoTc 


sev/r+Fv^r+e 


J (0, -B) - J (e, -B) + J (^, s) - J (0, B) 


(26) 


where the f un ction J is dehned as 

J (x, y) — \/+ x'^ — 2/ In ( 2 ^/A'^ + + 2a:j 

— a/A^ + In 


2^/A^ + y^VA^ + x^ - 2xy + 2A^ 


X + y 


(27) 


In these last two eqnations, B^ = 1 + and A‘^ = it ~ the center of the galaxy, 

the mean intensity is given by 

J, (0) = -^0L= (28) 




and at the onter edge of the galaxy its valne is 


The local heating (ergs s“^) of a spherical dnst grain of radins a by ambient starlight 

'■y(^ 

Here c is the speed of light, Qi, is the dnst absorption efficiency (kindly provided by Dwek 
1993), and the integration is over optical and nltraviolet freqnencies. 

To evalnate i7*, we mnst determine the coefficient jo- Assnming that the stellar 
mass to light ratio is spatially constant, the ratio of the stellar mass density to the stellar 
emissivity in a given freqnency interval is also constant by assnmption (21). Therefore the 


cQjjTiordi'. 


(30) 
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ratio of p*o in eq. (1) to jo in eq. (21), when integrated over the B-band, is hxed by the 
stellar mass to B-band light ratio of eq. (A2), 


joB _ f 

p*o V J 


(31) 


The determination of JoB allows the total heating dne to starlight in the B-band to be 
compnted. However, we reqnire the heating dne to all of the starlight, not jnst the light in 
the B-band. This is obtained by mnltiplying the heating dne to starlight in the B-band 
with a bolometric correction factor given by 


J„V. optical F W QA<iA 

/b band ^ 


(32) 


where Q\{a) is again the dnst absorption efficiency at wavelength A. The npper integral 
extends over all optical and nltraviolet wavelengths, the lower integral extends only over 
the B-band, and F (A) is the spectrnm of the stellar light in ellipticals. We assnme that 
F (A) does not depend on the radins ^ and we nse the analytic £t 


' log (A) - 3.215 
-1.509 log (A) +4.314 
- 0.6 


log F{X) = { 


V 


4.078 log (A) - 14.17 

0.86 

-0.8812 log (A) +4.341 
-2.837 log (A) + 12.58 


if A < lOOOA 
if lOOOA < A < 1800A 
if 1800A < a < 2123A 
if 2123A < A < 4842A , 
if 4842A < A < 8913A 
if 8913A < A < 16255A 
if A > 16255A 


(33) 


where A is given in A and F is in arbitrary nnits; the absolnte normalization of F (X) 
is given by eq. (31). Eqnation (33) gives a good fit (< 15% in the optical part) to the 
spectrnm of stellar light in elliptical galaxies as determined by Brnznal (1985), Renzini & 
Bnzzoni (1986), and Fergnson & Davidsen (1993). 


5.2 Heating by Hot Electrons 

While dnst grains are spnttered by collisions with thermal protons and helinm nnclei, 
collisional heating of grains is dominated by more freqnent collisions with eqnally ener¬ 
getic bnt more rapidly moving electrons in the hot interstellar gas. At plasma temperatnres 
T > lO^K the electric charge of the grains has a negligible inflnence on these collision rates 
(Draine & Salpeter 1979). For larger grains the energy transferred from the colliding elec¬ 
tron to the dnst grain is small relative to the total thermal energy of the grain. Heating 
therefore occnrs in a smooth fashion, reqniring many electron impacts. For small grains, 
however, a single electron collision can introdnce enongh energy to significantly raise the 
temperatnre of the grain. For small grains, electron collisions are snfhciently infreqnent 
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that substantial cooling by radiative losses may occur between successive impacts. Al¬ 
though small grains are stochastically heated (Dwek 1986), the time-averaged heating rate 
due to hot electrons is given by eq. (3) of Dwek (1986): 


) [ Ef{E)viE)C{E)dE (34) 

V^p/ Jo 

where E and v (E) are the energy and velocity of the impinging electron, and p is the gas 
density. The energy distribution of the electrons is Maxwellian, / [E) = {kT) 

i?^/^exp (—i?//cT), where T is the gas temperature and k is Boltzmann’s constant. The 
function C,{E) gives the fraction of the energy of an incident electron that is absorbed by 
the grain (Dwek & Werner 1981), 


C{E) 


1 - {E,/Ef/^ 


2/3 


iiE<E^ 
ii E>E^^ 


(35) 


where if* (ergs) = 3.7 x 10 ® {a/. 


5.3 Heating by X-Rays 

Some fraction of the thermal X-rays emitted by the hot interstellar gas in ellipticals 
is absorbed by the dust. Since the energy contained in each X-ray photon is comparable 
to the thermal energy of a grain in equilibrium with the stellar radiation held. X-rays also 
cause stochastic heating of the grain. Stochastic heating can also occur for the smaller 
grains when an ultraviolet stellar photon is absorbed (see e.g., Draine & Anderson 1985). 
Although stochastic temperature excursions are important in considering the predicted 
emission spectrum from the grains, we restrict our attention here to the average heating 
rate due to X-ray absorption. 

The X-ray emissivity of the hot gas is given by 


jx (^, i') 


_ rup _ 


Ax (T,z/) 


(36) 


where Ax is the atomic emissivity of the hot gas (ergs cm^ s“^ Hz“^) and T is the gas 
temperature. Since p is given by eq. (4), jx can be represented by eq. (21) with the 
identihcation jo = {po/mpYKx (T, z/). The heating rate can then be computed using eq. 
(30) integrating over X-ray frequencies. We assume = 1 so that we actually compute 
an upper bound to the X-ray heating. The integral in eq. (30) reduces to an integral over 
only jo since the entire frequency dependence of J^, is contained in jo oc Ax (T, z/). 


5.4 Comparison of Heating Rates 

The heating rates per grain for each of the three processes discussed above are plotted 
as a function of grain size in Figure 4 for the specihc case of graphite grains at a radius 
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of 0.4 in galaxy model a. Heating due to starlight is orders of magnitude greater than 
that due to either electronic collisions or absorption of thermal X-rays. The dominance of 
heating by starlight absorption holds at ah radii for ah galaxy models in Table 1 and for 
both graphite and silicate grains. Heating rate variations with grain size for other radii and 
galactic models are qualitatively similar to those shown in Figure 4 and are not presented. 

These results are consistent with previous hndings based on less quantitative treat¬ 
ments (de Jong 1986, Jura 1986) and can be expected from elementary considerations. 
First, for even the brightest ellipticals, the optical luminosity (B-band) is stronger than 
X-ray luminosity by a factor of ~ 100 (see eq. [6]) and this factor carries over to grain 
heating. Further, unless the energy lost from the hot interstellar gas through electronic 
heating of grains is vastly larger than that lost through plasma radiation, grain heating by 
electronic collisions must also be small (the various channels of energy loss for the hot gas 
are considered in §7). 

Similar energetic considerations for cD galaxies or for elliptical galaxies at the centers 
of clusters of galaxies would be very interesting since it has been suggested that the princi¬ 
ple grain heating mechanism is electronic collisions rather than optical heating (Bregman, 
McNamara, & O’Connell 1990; de Jong et al. 1990). Based on their far larger ratios of 
Lx to Lb (e.g. Edge & Stewart 1991), cD galaxies have a much greater amount of hot gas 
per star then do isolated elliptical galaxies; they also have much larger ratios of Lib to Lx 
(Edge & Stewart 1991). The implication is that higher ratios of gas to stellar density in 
cD galaxies favor plasma cooling by electron-grain interactions. This will be considered in 
a subsequent paper. 


6 . EVOLUTION OF THE GRAIN SIZE DISTRIBUTION 
The equation for evolution of grains in a hot Lowing gas with sources is 


dN {a,r,t) d 

dt da 


N (a, r, t) 


da 

dt 


-|- V ■ [V (a, r, t)u] = S (a, r, t). 


(37) 


In this equation, N(a, r, t) is the grain size distribution at galactic radius r and is given in 
units of cm“^ium“^, d = da/dt is the sputtering rate given by eq. (14) and S is the source 
function for dust. Assuming spherical symmetry and steady state conditions, this reduces 
to 


- dN dN . 

-hp— -F u— = b (a, r) 


da 


dr 


JV|2“ + A 

r dr 


(38) 


where u is the radial velocity and h is given by eq. (19). We recognize that the source of 
grains S' cx a* (t) is a strong function of time, but it will be apparent from our results below 
that the steady state approximation is sufhcient to illustrate the effects of grain evolution. 
The quantity S(a,r) in eq. (38) is therefore the source function for dust at the present 
time, tfiow' 

The differential equation (38) is quasi-linear and can be solved by the standard method 
of characteristics. The velocity of the gas u{r) has already been determined by the cooling 
how calculations of §3, but the source function S and appropriate boundary conditions need 


14 



to be specified. Since we only consider dnst ejected from evolving stars, S is simply given by 
the rate of dnst inpnt from these stars. Unfortnnately, there are presently few observational 
constraints on the grain content of red giant winds since dnst absorption featnres are 
difficnlt to observe in the spectra of these cool stars. What constraints exist snggest that 
there are relatively few small grains (a < 0.08//m, Seab & Snow 1989). However, in 
view of the absence of information on the grain size distribntion we assnme the grain size 
distribntion is similar to that in onr Galaxy: a power law with an exponent given by the 
MRN (Mathis, Rnmpl, & Nordsieck 1977) valne of g = 3.5, 

S{a,r) = anP* (a~^, (39) 

\pmpj 

where A is an as yet nndetermined normalizing coefficient. We adopt a minimnm grain 
size amin = 0.001/um and a maximnm grain size of ttmax = 0.3/um. 

In the standard MRN model, the valne of A is determined by extinction measnrements 
of the Galactic interstellar medinm. Since these data do not exist for onr case, we determine 
A based on the local metallicity of stars in elliptical galaxies. Gonsider the grain size 
distribntion of an evolving star: 


^ej{,Oj) Pej-^ejOj 


(40) 


where p^j is the gas density in the ejecta. The total mass density of grains in the ejecta is 
then 


Pdust 


TTCi pgrain ) Pej-^ej^ ^dCL 


(41) 


where Pgrain is the mass density of an individnal grain. We take Pgrain = 2.2 g cm“^ for 
graphite and Pgrain = 3.0 g cm“^ for silicate. The density of dnst is related to the gas 
density by the dnst to gas ratio Pg: 


Pdust — VgPej- 


(42) 


Then, from eq. (41): 


A, 


ej 


Vg 

4 

3 ^ Pgrain 



(43) 


The solar abnndance by mass of both silicon and carbon is abont 0.005. We therefore set 


yg = 0.005 3 Z (r, Lb) , 


(44) 


where 5 is the fraction of either Si or G in the ejecta contained in grains, and Z is the 
metallicity of the stars relative to solar. Based on optical line width measnrements (Davies, 
Sadler, & Peletier 1993; Schombert et al. 1993; Gonzales et al. 1994), the metallicity is 
known to vary with radins and B-band Inminosity of the galaxy approximately as 


Z{r, Lb) 


1 , L / ^0•3 

1 + (r/rc) 


(45) 
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where Tc is the galaxy core radius. We therefore hnd that 


^ = A, =( ^ ^ 

l^'^p X'^'^PgrainJ \y(lrnax ^min / 


(46) 


where i/g is given by eq. (44). In deriving S from (39) we have assumed that few grains are 
sputtered before the stellar ejecta merges into the hot ISM. We provide a brief argument 
in Appendix C to support this assumption. 

It is useful to compare our values for the normalization A with standard values associ¬ 
ated with the MRN distribution. For carbon grains in a typical galaxy of lO^^L©, we hnd 
A = 0.053 at the core radius Tc- This assumes that all the carbon in the ejecta is incorpo¬ 
rated into grains {Sc — !)• This is close to the standard MRN value of A = 0.069 (Draine 
& Anderson 1985) for graphite grains in the Galaxy. For silicate, again assuming dsi = 1, 
we get A = 0.039 which is similar to the MRN value A = 0.078. Although our adoption 
of the source function in eqs. (39 - 46) is approximate, a detailed comparison of infrared 
observations with computed IRAS spectra will lead to an improved understanding of the 
grain abundance and size distribution in the ejecta of red giants and planetary nebulae. 

Regarding an appropriate boundary condition at large galactic radius for solving eq. 
(38), we note that a favorable aspect of our method is that the grain size distribution 
in the galaxy actually does not depend on the choice of the outer boundary condition 
although something has to be specihed for the numerical calculation. This is seen by 
referring to the timescales of Figure 3. In all cases, towards the outer boundary, the flow 
time far exceeds the sputtering time of the largest assumed grain. Indeed, the flow time 
tends to inhnity as the flow velocity decreases toward the outer boundary. The grains in 
the outer regions therefore sputter on the spot and the grain size distribution at small or 
intermediate galactic radii is insensitive to the distribution at larger radii. We may specify 
any grain size distribution at large radius to start the integration of the grain flow equation 
(eq. [38]). For consistency, we assume that the grain content at large radius is determined 
by the on the spot sputtering assumption of eq. (18) with the source function of eq. (39), 
i.e. 


N (a, r) 


{g - 1) hp 


A 


prUp 


, 1-9 


(ft ^ O'max 


1-9 


)• 


(47) 


Finally, we must specify that the number density of grains vanish for the biggest grains in 
the ISM ttmax, N {ajnax,r) = 0 at all r. 

Figure 5a shows grain size distributions for graphite grains determined from eq. (38) 
for model galaxy a. The grain size distributions are given at the core radius Tc (topmost 
curve), and at successively larger radii (lower curves). The hrst implication of this hgure 
is that the grain size distribution is very close to that of in situ sputtering (eq. [18]). For 
small grains the distribution is essentially a power law that is shallower than that of the 
stellar ejecta by unity. At large grain sizes, the density steepens and drops to zero at the 
largest assumed grain size. We anticipated this result wherever the sputtering times are 
shorter than the flow times (Figure 3). What is somewhat surprising, however, is that the 
assumption of in situ sputtering can be used to accurately determine the size distribution 
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even where the sputtering time for the largest grains is comparable to or exceeds the flow 
times to the core. This can be understood since the flow of grains from upstream into 
a given volume near the galactic core is signihcantly smaller than the amount of 

grains produced by stars within A'lrr'^dr - this is a result of the steep stellar density prohle 
p*(r). 

The grain size distributions for graphite in the smallest galaxy (model c) are shown 
in Figure 5b. In contrast to galaxy model a, the flow time to the core radius is comparable 
to the sputtering time for the largest grain over most of the galaxy. Again, the grain size 
distributions are well determined by assuming in situ sputtering. Similar results hold for 
galaxy model b which are not illustrated. Grain size distributions for silicate dust have 
identical shapes to those of graphite, although the normalizations are slightly different, so 
these results are not shown. A favorable consequence of the above results is that the grain 
size distribution in the ISM depends only on the ratio of the stellar density to the gas 
density (and on the rate of dust injection into the ISM). These quantities are reasonably 
well known. The grain size distribution does not depend on the more poorly known gas 
velocities in the cooling flow. 

In order to gain further insight, we have varied ttmin and ttmax while keeping the 
slope of the grain distribution g hxed. For example, if instead of Umax = 0.3/um we set 
(imax = 1.0/Um in the stellar ejecta, the grain size distribution is modihed in the ma nn er 
shown in Figure 5c. Over the common range of grain sizes, the case with Umax — 1.0/um 
gives about half the grain density relative to the case with Umax = 0.3/um. Since the total 
mass of grains ejected by stars is hxed, extending the range over which grain masses are 
distributed requires a decrease in the density per unit grain size. That is, by increasing 
dmax we decrease the normalization A by a factor of ~ 0.53. 

One hnal case we consider is motivated by the observational evidence that there may 
be fewer small grains in stellar ejecta than in the Galactic ISM. For example, consider a 
power law grain size distribution S{a) oc with a^in = 0.08/um < a < ttmax = 1-0/um. 
In the cooling how ISM, however, we continue to follow the evolution of grains smaller 
than ttmin since these smaller grains are produced by sputtering. The resulting grain size 
distribution is shown in Figure 5c as the dashed line. Since the range of grain sizes in 
the source function S'(a) is restricted at small grain radius a, the density per unit grain 
size is higher than in the case where ttmin = 0.001/um and comparable to that of the hrst 
case considered (solid line of Figure 5c). For a < ttmin the grain size distribution in the 
ISM is essentially constant. Because there are no stellar sources of dust in this range, 
the only sources are the sputtered remnants of larger grains. Since the sputtering rate is 
independent of the grain size (eq. [14]), all grains get smaller at the same rate. Similar 
results are obtained for galaxy models b and c but are not illustrated. 

7. ENERGY LOSS FROM THE INTERSTELLAR GAS 

Given the grain size distributions of §6, we can compute the total rate of energy loss 
from the hot ISM by electron-grain collisions. This can be done at any given radius by 
simply convolving the heating rate per grain (eq. [34]) with the grain size distribution. 
This rate of energy loss from the plasma by electron-grain collisions can be compared with 
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(48) 


the local rate of energy loss due to thermal emission, 

2 

Atot (T) , 


dE 


X 


dt 


p(0 


rur 


where Atot (T) is the bolometric atomic X-ray emissivity. 

Cooling rates for the hot gas in the largest galaxy (model a) are shown as a function 
of galactic radius in Figure 6. In all cases, the energy loss rate from the hot plasma due to 
grain heating is signihcantly smaller than that of the thermal emission. This result is very 
different from that of Dwek (1987) who found that dust heating accounted for as much 
two orders of magnitude more energy loss from hot, dusty plasmas than thermal emission. 
Dwek then applied these results to supernova remnants. The difference between our results 
and those of Dwek’s is our inclusion of sputtering losses. To estimate the magnitude of the 
reduction of the dust-to-gas ratio due to sputtering, assume that all grains have the same 
radius a and mass uig and that they are created and destroyed at the same rate. In this 
case the number density of grains N is related to the dust-to-gas ratio in the stellar ejecta 
i/g by Nuig/tsp = cxnP*yg, where tgp ~ arup/hp is the local sputtering time in plasma 
of density p (eqs. [14] and [15]). When sputtering is included the dust-to-gas ratio is 
Vgs = 'nT'gN/p ~ VgtspCinP*/P- It follows that sputtering reduces the dust-to-gas ratio by 
a factor Pgs/Vg ~ igpOinP^/p ~ {ampan/h)p^,/p‘^. Consider for example stellar and gas 
densities near the core radius of galaxy model a, p* ~ 3 x 10“^^ g cm“^ and p/rrip ~ 0.1 
cm“^. At this location the dust-to-gas ratio is reduced by Pgs/Vg ~ 3 x 10“^ for grains of 
size a = 0.01 pm; this is essentially the factor that accou n ts for the difference in our plasma 
cooling rates and those of Dwek (1987) who assumed an “extended” MRN distribution with 
an unsputtered dust-to-gas mass ratio. By allowing the MRN distribution in the stellar 
ejecta to be sputtered, we hud that the plasma energy loss to grain heating is less, not 
very much greater, than radiative losses. We also hud that the space density of the largest 
grains is preferentially decreased; these are the grains most responsible for energy loss from 
the hot gas (see below). Self-consistent grain sputtering therefore plays a decisive role in 
correctly assessing the signihcance of dust emission for the hot ISM in galactic cooling 
flows. 

It is evident from Figure 6 that the heating of graphite grains (short dashed-dotted 
line) accounts for a little more than half of the total plasma cooling rate whereas the 
less numerous silicate grains (long dashed-dotted line) contribute somewhat less than half 
of the total. We also hud that the various assumptions made concerning the grain size 
distribution of the stellar ejecta do not result in large differences in the rate of heat loss 
from the hot ISM. This is easy to understand by considering the dependence of the plasma 
cooling rate on the assumed grain distribution of the stellar ejecta. The heating rate for 
the dust is 


dE. 


dust 


dt 


N (a, r) H^da oc 


4—g 
Ojmax d 


4-g 



^a^da 


Constant. (49) 


The term in parenthesis is due to the normalization of the grain size distribution (eq. 
[46]), the hrst factor in the second integrand is roughly the power law of the grain size 
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distribution in the hot ISM, and the second factor in the second integral is approximately 
proportional to the cross-section of the grain for electron-grain collisions. The plasma 
cooling rate due to electron-grain collisions is therefore insensitive to the assumed value 
of the slope g or the upper and lower cutoffs in the grain size distribution in the stellar 
ejecta. Eq. (49) shows that the plasma cooling rate depends only on the total mass of 
grains ejected by stars into the ISM and the ratio of the stellar mass to gas mass. Since 
these are well constrained by observations, the rate of plasma cooling by this process is 
well determined regardless of specihc assumptions made about the grain size distribution 
in the stellar ejecta. 

According to the results plotted in Figure 6, the ratio of energy loss from the hot 
ISM by dust heating to that by X-ray emission is relatively independent of galactic radius 
r. The heating of individual grains by hot ISM electrons is proportional to the local gas 
density, o: p (eq. [34]). But the local grain density depends on the ratio of stellar to gas 
density and the stellar metallicity, N{a, r) oc Zp^j p\ therefore NH^. oc Zp^. By comparison 
the rate that the hot ISM loses energy by thermal radiation dEx /(it oc oc p* (eqs. [1], [4], 
and [48]). Therefore the ratio of thermal losses by grain heating to optically thin emission 
varies slowly with galactic radius, NHe/{dEx/dt) (x Z(r) oc which accounts for the 

small rate of divergence between the two rates plotted in Figure 6. For the largest galaxy, 
for example, (r^/rc)®'^ ~ 159^'^ = 4.6 is the ratio of the two heating rates at the core 
radius and at the outer radius, again consistent with Figure 6. 

The loss of thermal energy from the hot ISM by grain heating in the two smaller 
galaxies b and c is even less important relative to radiative losses than in model a. This 
can be understood because NH^/{dEx/dt) oc Z oc Lb (eq. [45]). The larger ratio of 
sputtering to flow times in smaller galaxies (Figure 3) plays a very minor role in this ratio; 
the relative contribution of advected dust to the total dust content at any radius is small 
because of the steep radial dependence of p* (r). 

8 . CONCLUSIONS 

We have investigated the evolution and energetics of dust embedded in the hot ISM 
of isolated elliptical galaxies. For this purpose we constructed a series of single-parameter 
galaxy models assuming King type prohles for the stellar density which are consistent with 
the fundamental plane relation, the mass to light ratio, core radius relation, and a virial 
condition for isothermal cores. We consider three representative model ellipticals with Lb 
spanning a range from IQiOLo to IQiiLo which are expected to contain galactic 
cooling flows. 

The galaxy models were used to construct approximate isothermal - but fully time- 
dependent - cooling flow models for the hot ISM. The radial gas density variation is known 
if the X-ray surface brightness distributions of all ellipticals resemble those resolved by 
Einstein. With this assumption we can solve for the velocity held using the continuity 
equation alone. These “continuity” models compare well with detailed hydro dynamical 
cooling how calculations. The isothermal temperature of the gas is taken to be the same 
as the mean stellar temperature of the galaxy. This latter temperature is determined by 
solving the equation of stellar hydrodynamics in the presence of a massive dark halo. These 
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temperatures agree reasonably well with observed temperatures for galaxies of comparable 
luminosity (see Figure 1). A more detailed agreement of predicted and observed tempera¬ 
tures is not required for the purposes of this paper because the results depend only weakly 
on the assumed temperature of the ISM. 

Having specihed the plasma enviro nm ent within the model galaxies, the radiation 
enviro nm ent from the near IR through the UV is taken from observations. With this in¬ 
formation we can compute the heating rate of interstellar grains due to absorption of stellar 
light, impacts with energetic electrons from the hot ISM, and absorption of thermal X-rays 
from the hot ISM. In spite of the efficiency of sputtering in grain destruction, the heating 
of grains in all model galaxies is dominated by absorption of ambient starlight by over an 
order of magnitude. This result applies for both graphite and silicate grains. Although 
this result was anticipated by de Jong (1986) and Jura (1986), our detailed calculation 
of the heating rates is relevant to the well constrained radiation-plasma enviro nm ent of 
elliptical galaxies. In future work, we will consider the same issue for cD galaxies and giant 
ellipticals at the centers of clusters of galaxies where both the plasma environment and IR 
emission characteristics appear to be distinct from isolated ellipticals (see §5.4). 

We have computed the grain size distribution in three representative galactic cooling 
flows including the effects of grain sputtering, advection by the cooling flow, and local 
stellar sources of dust. In this paper we have not considered external sources of grains 
- for example by recent mergers of gas-rich galaxies with ellipticals - or dust that may 
form within gas deposited in the centers of the cooling flows. In this sense our estimate of 
the amount of dust present in ellipticals is a lower limit. We assumed that the grain size 
distribution introduced into the cooling flow at every radius from evolving stars is given 
by a power law having the MRN index of g = 3.5, normalized according to the observed 
local metallicity of the parent stars. No signihcant loss of grains is expected during the 
brief transition from red giant winds or planetary nebulae into the ISM. The ratio of grain 
sputtering time to radial flow time varies within each galactic cooling flow and among 
galaxies of different total mass. When the sputtering time is very short we hnd that the 
grain size distribution for small grains will be a shallower power law with index 1 — g, but 
we hnd that this same solution is still an excellent approximation even when the sputtering 
time is so long that the grains can move inward during their lifetimes across an appreciable 
part of the galaxy. This last more surprising result can be understood from the very steep 
radial gradient of the stellar sources of grains; grains advected from larger radii in the 
galaxy are far less numerous compared to those produced locally at every radius in the 
cooling how. This latter property implies that the grain size distribution does not depend 
on the huid velocities in the cooling how. Rather, the distribution depends only on the 
better known ratio of the stellar density to the gas density. 

The heating of grains by absorption of hot thermal electrons is a possible channel of 
energy loss for the hot gas in cooling hows. Dwek (1987) suggested that this mode of energy 
loss may far exceed that due to emission of thermal X-rays in hot plasma environments 
similar to that of the ISM of ellipticals. However, we hnd that this channel for energy 
loss is in fact much smaller than that due to thermal emission in galactic cooling hows. 
While Dwek’s result was based on assuming that the grains (with an MNR distribution) 
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are unsputtered while they reside in the hot plasma, in our model galaxies the dust-to-gas 
ratio and the distribution of grain sizes are considerably modihed by sputtering. The total 
plasma cooling by a combination of graphite and silicate grains is a factor > 5 below that 
of optically thin radiative losses in all galaxies we consider. Our results for the magnitude 
of plasma cooling by thermal interactions with grains are very insensitive to the assumed 
upper and lower grain size cutoffs as well as the slope g of the grain distribution in the 
stellar ejecta. 

In this paper we have surveyed the relevant physics of distributed grains in elliptical 
galaxies. A second paper, currently in preparation, will describe in detail the temperatures 
of grains at every radius in elliptical galaxy cooling flows and the far infrared emission into 
each of the IRAS bands. We will determine the radial distribution of IR emission for com¬ 
parison with future, high resolution observations. We are also interested in estimating the 
IRAS luminosities from more centrally concentrated dusty patches and lanes in ellipticals 
for comparison with the more distributed dust. 

We are grateful to Eh Dwek for providing us with grain absorption coefficients, to Scott 
Trager and Sandra Faber for insights on observations of ellipticals, and to John Raymond 
for use of his code. JCT was supported by an NRC associateship and WGM thankfully 
acknowledges support from CalSpace and a Faculty Research Grant from UGSG. 
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APPENDIX A 

SELF-CONSISTENT GALAXY MODELS 

The fundamental plane for elliptical galaxies 

Lb = 3.6 X 10^ (^1) 

relates the B-band luminosity Lb (in units of Lqb), the effective radius Re (in parsecs) 

and the central stellar velocity dispersion a*o (in km s“^); h = H/lOO Mpc km“^ s“^ is 
the Hubble parameter (Dressier et al. 1987; Djorgovski & Davis 1987; Ciotti et al. 1991). 
The stellar luminosity Lb can be directly related to the total stellar mass (in units of 
Mq) by the stellar mass to light ratio: 

= 2.98 X 10“^ (A2) 

based on the recent study of van der Marel (1991). 

The effective radius Re that contains half of the stellar mass or light in projection 
can be found from eq. (3) in terms of Ve and rt- The resulting implicit equation for the 
effective radius Ae = Re/fc is approximated very well by 

Ae = 0.857 (A3) 

where = Tt/re- In addition, the core radius rc, in parsecs, increases with galactic 
luminosity according to 

Te = 6.84 X 10“^^ L)j20^1.4 

(Lauer 1989). This relation is consistent with more recent galactic distances based on the 
Dn — cr* method (Donnelly, Faber, & O’Connell 1990). 

Finally, the core radius, central density, and central velocity dispersion are related via 
the standard core relation for isothermal spheres (Binney & Tremaine 1987), 
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(A5) 


By elimination among eqs. (Al)-(A5), all quantities of interest can be expressed in terms 
of a single parameter, most conveniently chosen to be the total radius For a given 
the optical luminosity (in Lqb) is given by 

Lb = 8.92 X lO^' (^6) 


where Pmass is given by eq. (2) of the text. The central stellar density is then given by 
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The corresponding core radins and total stellar mass are fonnd from eqs. (A4) and (A2), 
respectively. Table 1 lists the parameters of a representative sample of the galaxy models 
generated by the above procednre. We assnme h = 0.75 thronghont. 

We now estimate the mean stellar temperatnre for onr galaxy models. Gas ejected 
from stars is heated to the virial temperatnre by dissipation. The gas can be fnrther heated 
by Type la snpernovae if the snpernova rate is snfhciently large. The observed snpernova 
rate is nsnally given in nnits of SNU = 1 snpernova every 100 years in a galaxy of = 10^° 
Lqb- For a typical large elliptical galaxy snpernova heating dominates heating by stellar 
motions if SNU > 0.1. However, recent observations of NGC 1399 and NGG 4472 by 
Serlemitsos et al. (1993) and Forman et al. (1993) indicate iron abnndances in the cooling 
flow gas that are near solar, similar to that of the stellar ejecta. The implied valne of SNU 
is less than 0.05 when compared with the models of Loewenstein and Mathews (1991). If 
snpernova heating is nnimportant, as seems likely, the gas temperatnre shonld be close to 
the mean stellar temperatnre (T*). We determine approximate valnes for (T*) by solving 
the eqnation of stellar hydrodynamics for the velocity dispersion (Mathews 1988) with 
plansible assnmptions for the dark halo mass distribntion and stellar velocity ellipsoids. 
The stellar temperatnre is then determined from the velocity dispersion by 


:(0 


limp 


3 GMc* 

2 Tc 




(A8) 


where iimp is the mean mass per particle (// is taken to be 0.5 and mp is the proton mass). 
Me* = Aur'^p^o/^ is the galactic core mass, and r is a dimensionless temperatnre. 

We assnme that dark halos in ellipticals are psendo-isothermal 




(A9) 


with a nniversal ratio of halo to stellar core radii /3 = Tchl^c = 10 and having the same 
onter radins as the stellar distribntion. The total integrated masses of the dark halos 
is assnmed to be larger than M*t by a factor w = 10. This implies valnes of pho listed in 
Table 1. For the stellar velocity ellipsoid we adopt a radial dependence given by 


m = 1 - 


■^tv 



(AW) 


(Mathews 1988) where and are the transverse and radial velocity dispersions respec¬ 
tively. 

The stellar temperatnres weighted by stellar mass and listed in Table 1 are compnted 
from the temperatnre distribntion T* (^) nsing 


(T*) = T*oT 
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We find that the average temperatures computed this way imply the relation (T*) cx 
or equivalently, Lb o: Since this is close to the observed Faber-Jackson law, 

we conclude that the assumed properties for the massive halo and velocity ellipsoid are 
reasonable. In addition, recently measured temperatures from ASCA (Serlemitsos, et al. 
1993; Mushotzky et al. 1995) agree well with eq. (All) and our relation (T*) oc 
(Figure 1). 


APPENDIX B 
THE L„,Lb relation 

Canizares, Fabbiano & Trinchieri (1987) showed that the X-ray luminosity in the 
EinsteinIPC band (AE ~ 0.5 - 4.5 keV) scales with optical luminosity as Lx oc Lb 
Later Donnelly, Faber, & O’Connell (1990) found Lx oc Lb using a smaller sample 

but with better galactic distances {Dn — method). Sarazin (1987) showed that Lx oc 
provided the cooling flow gas is heated by compression in the gravitational held; 
from this and the Faber-Jackson relation it follows that Lx oc L^^. In the following we 
present an alternate derivation of this result. 

Suppose that the hot interstellar gas in all bright ellipticals has a homologously similar 
distribution. If n is the mean plasma density then the rate that energy is emitted into the 
Einstein band per unit volume is proportional to u^Aae (Aas is the atomic emissivity 
in the Einstein band). The specihc emissivity is oc hAae /where iinip is the mean 
mass per particle. The total observed X-ray luminosity should be Lx oc nA^EMg/(imp 
where Mg is the total mass of hot gas. If (a*) is the specihc mass loss rate from an 
evolving population of old stars (Mathews 1989) averaged over a cooling how time (several 
Gyr), then by continuity of gas how Mg/tcool ~ {a^)M^f The cooling time is given by 
tcooi ~ ^kT/2Atotn where Atot is the bolometric cooling function. 

Combining these results, we expect 

Lx oc {a^){AAE/Atot)M^tT oc (a*)(AAE/Atot)M*iL|^’^, (51) 

where we have used the Faber-Jackson law Lb oc a* and the virial condition that (7* oc 
(T*) ^ T. 

If the stellar population in all ellipticals is dominated by equally old stellar populations, 
a^(tn) should not vary with Lb- However, eq. (Bl) involves (a*(t)) which is averaged over 
the time for gas to how through the central regions (^ < Ae) where most of the X-rays 
are produced. Typically this time is very short (< 10® yrs; see Fig. 2) so we conclude that 
(a*(t)) ~ a*(tn) does not vary appreciably with as long as all ellipticals have about 
the same age. For plasmas with solar abundance, Aae/A tot ~ 0.6 and is independent of T 
provided T > 3 x 10® K; Aae/A tot should not vary with Lb even for the least luminous 
galaxy in Table 1. However, the stellar mass to light ratio varies as M^t/LB oc where 
u ^ 0.2 (Dressier et al. 1987; Djorgovski & Davies 1987; Bender, Burstein, & Faber 1992; 
Renzini & Ciotti 1993) or n ^ 0.35 ± 0.05 (van den Marel 1991). Using this latter value 
and ?7 = 4 we hud Lx oc very close to the observed relation. 
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Finally we note from eq. (Bl) that Lj;{t) oc (a*(t)) oc which acconnts for the 

slow secnlar decrease in X-ray Inminosity dnring the cooling flow evolntion of ellipticals 
(Loewenstein & Mathews 1991; Ciotti et ah 1991). 

APPENDIX C 

DUST GRAINS IN STELLAR EJECTA ENTERING THE COOLING FLOW 


The complex evolntion of gas ejected from red giants as winds or planetary nebnlae 
as it interacts with ambient hot cooling flows has been discnssed by Mathews (1990). A 
typical evolntionary seqnence is: attainment of hydrostatic eqnilibrinm with the hot gas, 
disrnption by Rayleigh-Taylor instabilities into many small clonds that extend back along 
the stellar orbit as the orbital motion inherited from the parent star is dissipated, and 
condnctive evaporation into the hot ICM after coming to rest. Gas ejected from stars 
interacts only with the hot gas (not with other stellar ejecta) and can be kept ionized 
by galactic UV radiation prior to condnctive evaporation. However, gas ejected from 
stars mnst be condnctively heated to cooling flow temperatnres T ~ lO^K in less than 
~ 10® yrs or the collective emission line Inminosity wonld exceed that observed in elliptical 
galaxies. Condnction of thermal energy from the hot ISM into the stellar ejecta can occnr 
by magnetic reconnection as soon as the debris of the stellar ejecta slows below the local 
Alfven velocity. 

In this paper we assnme that dnst grains snrvive this transient evolntion from stellar 
ejecta into the cooling flow gas. The most convincing argnment for this is to adopt the 
opposite point of view, a worst case scenario in which the grains are destroyed before 
enconntering the hot ISM. Snppose, for example, that gas ejected from a star remains 
intact (i.e. no instabilities!) and that the density and temperatnre in the stellar ejecta are 
snfficient to allow spattering destrnction of grains before the gas density drops to ambient 
ISM valnes. Since spattering is very slow nnless T > I0®K, we snppose that the stellar 
ejecta, while still rather dense, is heated by condnction to the ambient temperatnre Ta 
of the local cooling flow. In this (very nnlikely) circnmstance the grains conld spatter 
away before entering the hot ISM. We now consider the self-consistency of this worst case 
evolntion. 

For plasma temperatnres T > I0®K the spattering time tsp for grains of size a = I0“® 
cm depends mostly on the plasma density in the clond tt - c , tgp = a/a ^ 1 x 10 ®a_ 5 n“^ 
yrs, where a _5 is the grain size expressed in nnits of 10“®cm (tt-c is in nnits of cm“^). A 
spherical clond of radins rc with the mass of an ejected stellar envelope, M^j = 0.2 M© and 
temperatnre Tc will expand into the low pressnre ISM in approximately the sonnd crossing 
time, tsc — Tc/cs ~ 2.5 x lO^Tcr yrs (Tcy is the temperatnre expressed in nnits 

of lO^K) where the clond is assnmed to have been heated to temperatnres comparable 
with the local ISM. In order for the grains to spatter away before entering the local ISM, 
we reqnire tgc > Gp or ric > 246a_5 cm“^. Since the internal pressnre in snch 

a clond, ricTc ~ 10^ cm“^K, far exceeds that in typical cooling flows, UaTa ~ 10^ cm“^K, 
the spherical clond of stellar ejecta will indeed expand. 
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However, the thermal energy necessary to heat the stellar ejecta to temperatnre 
mnst come from the local ISM. The radins of a region in the local ISM that contains the 
same thermal energy as the clond is = rdncT ^/~ 3.8 x ^^^{naTa)^ 

cm, where {naTa)^ is {riaTa) given in nnits of 10^cm“^K. Finally, thermal heat mnst be 
condnctively transported into the stellar ejecta before the clond can expand nnder its own 
pressnre. Therefore in order to destroy the dnst before it enters the local ISM we reqnire 
that heat is condncted into the spherical clond with a characteristic velocity Ufi > Tajisc ~ 
4.8 X lOV 5/6(n„T„)4 -1/3 km s”! > 3 x lO^a.s “^/^ km s”!. 

Since this lower limit exceeds the velocity of thermal electrons in a gas at ~ lO^’K, 
fe ~ 2 X lO^Ty 1/^ km s”!, we conclnde that this hypothetical, maximally pessimistic 
evolntionary model is impossible. Dnst grains do in fact largely snrvive the transition 
from stellar ejecta into the local cooling flow ISM as we have assnmed here. 
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TABLE 1 


Galaxy Model Parameters^ 


Model Lb{10^^Lq) 

rc(kpc) 

n(kpc) 

P*o 

0 ^ 

Ho 

„ b 

Pho 

(T*)(107K) 

a 10.6 

0.774 

123. 

3.37 

2.42 X 10“^ 

1.12 X 10-2 

0.922 

b 3.31 

0.192 

76.4 

38.6 

5.06 X 10“^ 

5.74 X 10-2 

0.468 

c 0.976 

4.44 X 10-2 

44.4 

518. 

1.47 X 10-3 

3.48 X 10-1 

0.252 


“ Parameters for the galaxy models given in the hrst colnmn are listed. From the 
second to the eighth colnmn, respectively, these are: the B-band Inminosity, the core 
radins, the total radins, the central stellar density, the central gas density, the central dark 
matter halo density, and the mean temperatnre of the stellar distribntion. 

^ All densities are given in nnits of 10“^^ g cm“^. 
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FIGURE CAPTIONS 


Figure 1: The (LB,(T*))-relation for the model galaxies of Appendix A is shown. 
The hlled circles (connected by the solid line) give resnlts for some of onr models. The 
open circles give the observed T and Lb for NGC 4636, NGC 4406, and NGC 4472. 

Figure 2: The solid line gives the inward velocity of the gas in the cooling flow of 
galaxy a (see Table 1) plotted as a fnnction of the radins normalized to the onter radins. 
The dashed and dashed-dotted lines give the velocities for galaxies b and c, respectively. 

Figure 3: The sputtering time (dehned by eq. [15] of the text) and flow time (eq. 
[16]) are shown by dashed and solid lines, respectively. The heaviest lines give the results of 
galaxy a, the second heaviest lines correspond to galaxy b, and the lightest lines correspond 
to galaxy c. 

Figure 4: The heating rate (plotted versus grain size) due to absorption of ambient 
starlight is given as the solid line, the heating due to collisions with energetic electrons 
is given by the dashed-dotted line, and an upper bound for heating due to absorption of 
thermal X-rays is given by the dashed line. The indicated rates are computed at a radius 
of 0.4 times the outer radius of galaxy a (see Table 1). 

Figure 5a: The grain size distribution (graphite) at various galactic radii are shown 
for model galaxy a. The grain size distribution extends from amin = 0.001/um to amax = 
0.3/um. The topmost curve gives the number density of grains per unit grain size at the 
core radius The ten lower curves starting from the top give the grain size distributions 
at 0.1, 0.2,..., and 1.0 times the outer radius r^. 

Figure 5b: The grain size distributions (graphite) for model galaxy c are shown at 
various radii. The topmost curve applies at the core radius, and the ten successively lower 
curves correspond to 0.1, 0.2,..., and 1.0 times the outer radius. 

Figure 5c: The grain size distributions (graphite) at 0.5 times the outer radius for 
galaxy model a are shown. The solid line gives the case where we assume amin = 0.001/um, 
dmax — 0.3/um, and that the grain size distribution in the stellar ejecta is given by eq. (39). 
The dashed dotted line corresponds to the same case as the solid line except we assume 
o^max — 1.0//m. The dashed line assumes the same power law in the stellar ejecta, but that 
Q'min = 0.08/um and Umax = l-Ofim. The grain size distribution in the ISM is, however, 
computed for smaller grains than amin- 

Figure 6: The rate of energy loss from the hot ISM due to heating of grains is plotted 
as a function of radius for galaxy model a. The light solid line gives the total heating rate 
assuming amin = 0.001//m and amax — 0.3//m. The contributions to the grain heating rate 
from graphite and silicate grains for this case are given by the short dashed-dotted and the 
long dashed-dotted lines, respectively. The short dashed line gives the total heating rate 
in the case where amin = 0.001//m and amax = 1.0//m. The long dashed line corresponds 
to taking amin = O.OSfxm and amax = 1-Ofxm. The total heating rates given by the light 
solid line, the short dashed line, and long dashed line are nearly identical. The rate of 
energy loss due to X-ray emission is given by the heavy solid line. 
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